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Abstract 

We revisit the phase diagram of the = 4 SU{Nc) super- Yang-Mills theory coupled to Nf 
fundamental "quarks" at strong coupling using the gauge-gravity correspondence. We show that 
in the plane of temperature v.s. baryon chemical potential there is a critical line of third order 
phase transition which ends at a tricritical point after which the transition becomes first order. 
Close to the critical line there is an intriguing logarithmic behavior, which cannot follow from a 
mean field type of analysis. We argue that on the string theory side the third order phase transition 
is driven by the condensation of worldsheet instantons and that this transition might become a 
smooth crossover at finite 't Hooft coupling. 
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I. INTRODUCTION 



Fascinating new dynamical phenomena can appear when the couphng(s) of a system 
becomes strong. Famihar examples include color confinement and dynamical chiral symme- 
try breaking in QCD and high Tc superconductivity in various condensed matter systems. 
Yet strongly coupled systems are hard to solve and intuitions gained from weakly coupled 
systems normally do not apply. 

The AdS/CFT correspondence [1] has provided important tools for studying many 
strongly coupled systems by relating them to classical gravity (or string) systems. A re- 
markable feature of the duality is that highly dynamical, strongly coupled phenomena in 
gauge theories can often be understood on the gravity side using simple intuitive geometric 
pictures. For example, confinement is refiected in the fact that a fundamental string, which 
represents an unbounded quark in the boundary gauge theory, has no place to end in the 
bulk geometry and thus cannot exist. New strongly coupled phenomena may be waiting to 
be discovered, and AdS/CFT maybe just the tool needed to discover it. 

In this paper we are interested in understanding the phase structure of a large N^. gauge 
theory coupled to a small number Nf fundamental quarks at strong coupling from gravity. 
More precisely, we consider A/" = 4 SYM theory with a gauge group SU (Ac) coupled to Nf 
(A/" = 2) hypermultiplets in the fundamental representation of SU{Nc), which in the limit 
of Nc ^ Nf and strong 't Hooft coupling can be described in terms of A^^ probe D7-branes 
in the AdSs x 5*5 geometry [2]^. Comparing to QCD, the system has the following distinct 
features: 

1. Nc > Nf, while in QCD A^, ~ Nf. 

2. There are both fermionic and bosonic "quarks", which are charged under a U{1)b 
baryon number symmetry. 

3. The beta-function for gauge coupling is zero to leading order in Nf/Nc- The scale of 
the system is set by the quark mass m^. 

Despite many important differences from QCD, the system appears to be rather interesting 
in its own right and provides a nice laboratory for studying strongly coupled quark-gluon 

^ See [3] for a review of the system. 
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systems under extreme conditions. At zero baryon density and finite temperature it has 
been used to model heavy quark mesons in QCD [4, 5]. At finite baryon density and zero 
temperatures some novel dynamical features were recently found in [6]. 

The phase diagram of the system at finite temperature T and baryon chemical potential^ 
fig can be worked out by studying possible configurations of D7-branes in the background 
geometry (which is a black hole in AdSs x S^) and has been studied by various authors 
in [7-14] ^. The results can be summarized as follows (see Fig. 1): 

1. There is a transition curve in the fig — T plane which intersects with the horizontal 
axis at fig = nig with rrig the bare quark mass, and with the vertical axis at some 
temperature T = T^. Except for a small region near the vertical axis, the transition 
curve is given by 

N = ^?^\T) (1) 

where rrig is the effective quark mass at finite temperature (it decreases with tem- 
perature).^ In particular, the transition is second order [11] along the horizontal 
axis (T = 0) and first order along the vertical axis {fig = 0) [7]. The transition 
along the vertical axis at T = has been interpreted as a dissociation transition for 
mesons [7, 19, 20]. 

2. The region inside the curve in the fig — T plane is described by a D7-brane embedding 
which lies entirely outside the black hole (Minkowski- type embedding), while the region 
outside the curve is described by a configuration in which part of the D7-branes falls 
into the black hole (black hole- type embedding). In terms of the boundary gauge 
theory, the baryon number density is zero inside the curve and becomes nonzero outside 
(except along the vertical axis fig where the baryon density is always zero.) 

In this paper we improve on the above description and show that the transition along the 
curve (1) is a continuous (3rd order) phase transition which connects to a first order line 

^ In this paper we will use the terms baryon charge density and quark charge density interchangeably 
without a factor between them. In other words we take the U{1)b charge of a quark (anti-quark) 

to be 1(-1). 

For related studies of the same system see [15-17] and for related studies of a confining theory see [18] 
^ For a precise definition of TOg"^'' see discussion around eqs. (12)-(14) in sec. II A. 
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FIG. 1: The phase diagram in the ^q — T plane. The dashed line indicates a continuous transition. 
The transition line lies exactly on the curve rriq until a critical temperature very close (but 
not equal) to the dissociation temperature T^. 

(near the vertical axis) through a tricritical point whose location we identify precisely, see 
Fig .1. 

While our conclusions of a continuous phase transition along (1) are consistent with the 
analysis of [10] who noted a phase transition along the mq line at zero density they are 
different from a later discussion in [12] where a first order transition was noted. In [12] the 
relation between charge density and temperature was studied numerically near the transition 
line and a very small discontinuity in density was noticed, due to a change in dominance 
between a Minkowski type embedding with zero density and a black hole embedding. It 
appears to us that the discontinuity likely has to do with the numerical accuracy of their 
calculation. We do not find a discontinuity in charge density along the transition line (see 
equation (2) below) ^. More importantly, we have identified a clear physical reason for 
the phase transition, which indicates that it should not be considered as an exchange of 
dominance between different embedding solutions, which was behind the reasoning of the 
conclusion in [12]. As observed in [21], for Minkowski type embedding at finite temperature, 
there are worldsheet instanton corrections to the leading order Dirac-Born-Infeld (DBI) ac- 



We have obtained our results both analytically and numerically with agreement. 
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tion for the D7-branes. While for a Minkowski type embedding, particle excitations on the 
D7-branes are not sensitive to the value of the chemical potential, the worldsheet instantons 
which correspond to semi-classical strings stretched between D7-branes and the black hole 
do. In particular, when the baryon chemical potential exceeds the value (1), the instanton 
actions exponentially dominate over the DBl action in the large 't Hooft coupling limit and 
induce an instability. The consequence of the instability is that instantons "condense" and 
generate a genuine neck between the D7-branes and the black hole. Thus beyond (1) the 
Minkowski-type embedding cannot exist and should be replaced by a black hole type em- 
bedding.^ The transition also has a simple interpretation in the gauge theory. As discussed 
in [21], the instantons can be interpreted in boundary gauge theory as thermal medium 
quarks. When fig exceeds the value of (1), the quarks have negative free energies and will 
condense and generate a finite charge density, although it is not clear whether it is fermionic 
quarks or bosonic quarks which are condensing. 

More explicitly, we find along the critical line (1) (approaching it from above), the baryon 
charge density and the chemical potential are related by 

fig - mf ) = -B{T)e log e + A{T)e + O(e') (2) 

where e is the quark charge density (normalized to be dimensionless) and A{T),B{T) are 
some functions of temperature. In the zero temperature limit, B(T) goes to zero and A{T) 
goes to a finite constant, and one recovers the second transition at T = (found in [11]), 
where various exponents are given by their mean field values. At any finite temperature, 
however, it is always the logarithmic term on RHS of (2) that dominates at small enough 
densities and as a result the transition becomes third order. The logarithmic behavior does 
not appear to have a mean field counterpart, and thus this is an example of continuous 
phase transition from gravity which does not obey the Landau-Ginsburg behavior. Such log 
terms however do appear to be a common feature of the renormalization group analysis of 
condensed matter systems at their upper critical dimension. 

There also exists a temperature Tc at which BiT^) = and beyond which B{T) < 0. At 

^ Note that this reasoning does not by itself imply that the transition should be continuous, it only indicates 
that the transition is not an exchange of dominance as is normally the case for a first order transition. 
Our explicit calculation shows that there is no jump in the baryon charge density across the transition 
line in the infinite 't Hooft coupling limit. 
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Tc the transition is again second order. For T > Tc, the transition becomes first order since 

(T) 

for a given /ig close to, but sma//er than rrig , now there are two black hole type embeddings 
with e 7^ in addition to the Minkowski type embedding. Connecting a continuous critical 
line and a first order transition line, the point (T = fJ-q = m!}[\T^) is thus a tricritical 
point. Again the critical behavior near the tricritical point is not the same as that of a 
Landau-Ginsburg type effective theory. 

While it is natural to connect the first transition near the tricritical point with that near 
the vertical axis at /Xg = 0, our approximation in (2) which applies to small densities does 
not extend all the way to /ig = 0. Thus that part of the phase diagram remains a conjecture 
at the moment. Also numerical work at finte density by other authors [12] suggest this 
should be the case, however a conclusive argument remains to be made. 

The phase diagram in Fig. 1 is for the \ = oo limit, where the quark charge density 
is identically zero inside the transition curve. At finite A, as pointed out in [21], for any 
/iq 7^ at any finite temperature, the baryon charge density is in fact nonzero, given by a 
dilute Boltzmann gas of quarks. This immediately raises the question whether the transition 
along (1) is an artifact of the infinite A limit and will be smoothed into a crossover at any 
finite A. To settle this question requires summing over the worldsheet instantons found 
in [21] and will be pursued in a separate publication. 

The plan of the paper is as follows. In section II we discuss the gravity description of the 
gauge theory system. In sec II A we review the general gravity setup for describing A/" = 4 
SYM theory coupled to Nf fundamental hypermultiplets at finite temperature. In sec. II B 
we review how to introduce a finite baryon chemical potential. In sec. II C we argue that 
the phase transition is driven by string worldsheet instantons. In sec III we give a detailed 
analytically derivation of equation (2) and verify it numerically. In sec. IV we discuss the 
thermodynamics of the system which can be gathered from equation (2). In the conclusion 
sec. V we discuss the connection of the critical point to the first order dissociation transition 
at /iq = 0. We also discuss what happens to the phase diagram at finite A. 
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II. GRAVITY DESCRIPTION OF THE GAUGE THEORY 



A. General set-up 

At finite temperature, J\f = 4 SYM theory with a gauge group SU{Nc) can be described 
by a string theory in the spacetime of a black hole in AdSs x 6*5, whose metric can be written 

as 

ds' = ^ {-fdf + dx') + + R'dnl (3) 
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where x = {xi, X2, 2:3) and 

dfll is the metric on a unit five-sphere 55. The string coupling gs and the curvature radius 
R (in units of a') are related to the Yang- Mills coupling qym and Nc by 

^T^9s = 9ym^ — = v^, a = gYAiNc . (5) 
a 

The temperature T of the YM theory is given by the Hawking temperature of the black 
hole, 

One can introduce "quarks" to the system by adding to A/" = 4 SYM theory A'^^ M = 2 
hypermultiplets in the fundamental representation of the gauge group. In the limit of large 
Nc with Nf finite this can be described in the dual string theory side by adding Nf D7- 
branes in the black hole geometry (3) and to leading order in Nf/N^, the backreaction of 
the D7-branes on the background geometry can be neglected. A fundamental "quark" in 
the YM theory can be described by an open string with one end on the D7-branes and the 
other end on the black hole. Strings corresponding to quarks and anti-quarks have opposite 
orientations. Open strings with both ends on the D7-branes can be considered as "bound 
states" of a quark and antiquark, thus describing meson-type excitations in the YM theory. 

In the limit 

Ac ^ 00, A ^ 00, Nf = finite (7) 

which corresponds to the limit (7^ ^ 0, ^ in the string theory side, the geometric 
embedding and the dynamics of the D7-branes can be described using the Dirac-Born-Infeld 
(DBI) action 

S = -NfT-r j d^i V-det + 27ra'F„,„) . (8) 



where gmn denotes the induced metric on the D7-brane and T7 = , — is the tension 
of a D7-brane. This will be the limit we work with for most of the paper. We will also 
consider in the discussion section what happens when one relaxes the X = 00 limit, where 
new interesting physics could emerge. 

To describe the embedding of the D7-branes in (3), it is more convenient to introduce a 
new radial coordinate u defined by 

' r^[u) = with no = ^ (9) 



in terms of which (3) can be written as 

ds^ = -!^f^u)dt' + ^-^dx' + ^{du' + u'dnl) 



where 



In (10), we have split the last term of the first line in terms of polar coordinates on x 
with dQ"^ denoting the metric on a unit three-sphere. The D7-branes can be chosen to lie 
along the directions C,"' = {t, x, 1^3, p) and using the symmetries of the problem the embedding 
in the two remaining transverse directions can be taken as 0(^°) = and ?/(^°) = y{p). y{p) 
can be found by solving the equation of motion obtained from the DBI action (8) with the 
boundary condition 

y{oo) = L = (2W )mq (12) 

where ruq can be interpreted as the (bare) mass of the "quarks". We denote the resulting 
embedding function as yo{p), which were first obtained numerically in [13]. At a small tem- 
perature, the brane lies entirely outside the black hole, as indicated schematically in Fig. 2. 
The brane is closest to the black hole at p = 0, where the three-sphere in last term of (10) 
shrinks to a point. Denoting 

yo{p = 0) = Lo (13) 

then the shortest open string connecting the D7-brane to the horizon has a mass in the YM 
theory 

(14) 



mf = ^/ dyVfq 



UO 



p=0 
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rriq can be interpreted as the effective mass of the "quarks" at temperature T. Note 

(T) 

that rriq decreases monotonically with T, since as we increase the temperature, the black 
hole becomes bigger and gravity attracts the brane more to the black hole. There exists a 
temperature Td (= 2.166mg/-\/A), after which the branes fall into the black hole (often called 
black hole embedding) through a first order phase transition [7, 20]. We will be interested 
in the temperature range smaller or of order in which regime we always have [(3 = 

(3m^P ~ O(yA) . (15) 




FIG. 2: An embedding of the D7 brane (green) in the AdS^ x black hole geometry which 
Ues entirely outside the black hole. The exact form of the embedding has been exaggerated to 
emphasize certain features. 



B. Finite baryon chemical potential 

The system possesses a U{1)b global baryon symmetry under which "quarks" transform 
nontrivially. The associated conserved current is dual to the f/(l) gauge field on 
the D7-branes. Turning on a chemical potential for the charge density in the boundary 
theory then corresponds to imposing the following boundary condition^ for 

Aq{P = Oo) = fiq . (16) 

The charged density Uq can be calculated from the bulk theory by 

Uq = U%p = oo) (17) 

''' Also, if the D7-brancs fall into the horizon there will also be a boundary condition that Ag = at the 
horizon. 



9 



where 11° is the canonical momentum conjugate to Aq (in terms of p-shcing) evaluated at the 
classical solution which satisfies the boundary conditions (12) and (16). Again by symmetry 
we will take Aq to depend only on p and the other components of the gauge field will be set 
to zero. 

The DBI action (8) can be written explicitly in terms of y{p), ^o(p) a-s 

S = -NfT, I d'^py/' [fq (1 + (y'r) - a''] ^ , (18) 

where prime denotes derivative with respect to p, and we have introduced 

a(p) = 27ra'Ao . (19) 

We will denote the integrand of (18) as C. Since (18) does not depend on a (only a') we 
have a conserved quantity e: 

e=^ = r (20) 

a" = fq{l+ y") ^4--2 ■ (21) 
From (17), we thus find the charge density in the boundary gauge theory is given by 

ng = NfT7{27ra'){2iT^)e (22) 

where the factor 27r^ comes from the volume of the three sphere. 

Since we will be expanding in terms of small density later, it will be convenient to scale 
coordinates and e so that they are dimensionless, i.e. 

u,y,P^ Lo{u,y,p), e-^Lle, a{p) ^ Loa{p) (23) 

where Lq is the location of the tip of the brane (13) before turning on a chemical potential^. 
From now on, u,y, p,e,a{p) are all dimensionless. After the scaling (11) become 



One can also choose to normalize them using L (12) which is more directly related to the field theory 
mass. But in our calculation below using Lq is slightly more convenient. Lq is fixed once the ratio rUq/T 
is given. Another alternative is to use the temperature but that will make the zero temperature limit 
more subtle. 
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The boundary conditions for y{p) and a{p) now are 



2'Ka' 



y{oo} 



a[oo] 



(25) 



and equation (14) becomes 



m 



(T) 



27ra' 



dy\/ fq 



.T 



(26) 



p=0 



To obtain the equation of motion for y it is convenient to perform a Legendre transfor- 
mation on (18) to express it in terms of e. The transformed action is 



H = ea + C = \fe 



(27) 



which leads to equation of motion 



y" 



3y' , 8 py'-y 



3y' 6 py'-y 



1 + y'^ p u^ri~^ — 1 + p^q^ 

When e = (i.e. zero density), equation (28) becomes 



p 1 + u'^rj 



4n-4 



(28) 



^ V ^ 8 py'-y 







(29) 



1 + 1/'^ p u^rj~^ — 1 

whose solution is given by yo(p) described earlier around equation (13).^ In Fig. 3 we plot 
how Lq/L and rj change with temperature for an embedding governed by yo{p). 



0.95 




T/r,, 




FIG. 3: Lq/L and rj for the Minkowski solution (e = 0) as a function of T/T^- 



As stated earlier we will only consider T <Td for which there is only one Minkowski embedding solution. 
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C. Phase transition driven by string worldsheet instantons 

For any value of fig (and T < To) , one can always have the following solution to equations 
of motion (21) and (28) with 

= /ig = const, e = 0, y{p) = ?/o(p) (30) 

for which (28) reduces to (29). This solution is somewhat peculiar, since it implies that the 
quark charge density Ug is zero even at a finite temperature and finite chemical potential 
fiq. This appears to contradict with field theory expectation that the density for quark and 
anti-quark should be given (at low density) by the Boltzmann distribution n± ~ e"^^*"' 
which gives a nonzero net charge density rig = n+ — n_ 7^ for a nonzero fig. There is in 
fact no contradiction, since (30) is the result in the X = 00 limit (supergravity limit) and in 
this limit due to (15), 

n± ~ rig ~ e~^, for fig < m^^-* . (31) 

The charge density is thus exponentially small in the large A limit and not visible to any 
order in the 1 / -\/A (or a') expansion. 

In [21], it was found for the embedding (30) there are non-perturbative open string world- 
sheet instanton corrections to the DBI action (8) which accounts for the exponentially small 
quark density. More explicitly, the instantons are given by open strings stretching between 
the tip of the D7-branes and the black hole horizon and winding around the Euclidean time 
direction, as indicated in the left plot of Fig. 4. They are classified a winding number n. 
From the spacetime point of view, these instantons generate tiny virtual necks which con- 
nect the tip of the branes to the black hole horizon. The total Euclidean action for the 
D7-branes, which gives the thermodynamic potential of the boundary gauge theory in the 
grand canonical ensemble, can be written as 

S = Sdbi + ^ /^n exp {-\n\P{m'^p - sgn[n]/ig)) H (32) 

n=±l,±2,--- 

where n sums over the instanton contributions, Dn arises from the worldsheet determinant 
for each instanton^^, and ■ ■ ■ denotes other perturbative 1/a/A (or a') corrections. These 
string worldsheet instantons have a simple interpretation in the boundary gauge theory as 

It also includes possible integrations over instanton moduli space. 
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representing thermal medium quarks. In particular, the n = ±1 terms in (32) are precisely 
what one expects of a dilute Boltzmann gas of quarks and anti-quarks.^^ 

Due to (15), when fig < m^"^^ the instanton contributions in (32) are exponentially small 
in compared with the DBl action Sdbi- B^^t foi' 

N > (33) 

the instanton sum will be exponentially large compared with the DBI action and the so- 
lution (30) can no longer be trusted. In particular, since the instanton contributions are 
dominating, we expect them to "condense" and create a genuine neck between the brane 
and the black hole. This implies in the range (33) a new solution to the DBI action with the 
branes going into the black hole should emerge, as indicated in the right plot of Fig. 4. Thus 

IT) 

in the infinite A limit we have a phase transition at /ig = niq where the Minkowski-type 
embedding (30) goes over to a black-hole-type embedding. The transition also has a simple 
interpretation from the gauge theory point of view; when /ig satisfies (33), the quarks have 
negative free energies and will thus condense and generate a finite charge density. 

The possibility for such a phase transition has been studied before in [9, 10, 12], where a 
single^^ black hole embedding solution was found for jiq > mi[^ and it was concluded in [12] 
that the transition was first order. One reason for the conclusion was that it appeared that 
the Minkowski embedding solution (30) was still valid for /ig > niq and thus the transition 
appeared to be a change of dominance between different solutions. As we discussed above, 
the Minkowski embedding solution should be replaced by a black hole embedding solution in 
the parameter region (33). Thus if there is only a single black hole embedding solution for 

(T) 

fiq > nig the transition is most likely to be continuous. Indeed we will show in section IV, 
the transition is third order for temperature not too high, but then becomes first order 
through a tricritical point}^ 

Higher n contributions should encode corrections due to Bose-Einstein or Fermi-Dirac statistics and other 
corrections due to interactions, 
at a temperature not too high 

Beyond the tricritical point it becomes possible to have multiple black hole embedding solutions for a 
given fig. See section IV. 
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FIG. 4: Left plot: For Minkowski type embedding (30), there exist worldsheet instantons which 
correspond to semi-classical strings stretching between the tip of the D7-branes to the black hole 
horizon and winding around the Euclidean time direction represented as the angle on the inset disk. 
The radial direction on the disk is the same as the radial direction in the y — p plane. Right plot: 

(T) 

for jiq > rriq instantons dominate over the DBI contributions and will condense to form a genuine 
neck between the branes and the horizon, i.e. one should have a black hole type embedding. 

D. Black hole embedding at finite density 

As discussed in last subsection, for fig > rriq we expect a black hole type embedding 
for D7-branes in the gravity side and a nonzero quark charge density on the gauge theory 
side. That the quark density for a black hole type embedding should be zero can be seen 
as follows. For a black hole type embedding, one should impose an additional boundary 
condition = for the gauge field at the horizon in order to ensure the regularity of the 
solution. From (16) this implies that Aq should evolve nontrivially from the boundary to 
the horizon, i.e. there is a nontrivial electric field in the radial direction on the branes, 
which from (17) in turn implies 7^ (and thus e 7^ 0). Thus we should now consider 
equation (28) with e 7^ 0. Conversely it has also been argued previously in [16] that for any 
e 7^ 0, only black hole type of embeddings are allowed since for a Minkowski type embedding 
it is not possible to have a non-singular distribution of a finite charge density on the branes^^. 

Let us now consider the boundary conditions for y{p) at the horizon for a black hole type 
embedding. We will denote the point that the brane enters the horizon and Uc = y{pc)- 
Note that 

yl + pI = ul (34) 
See also [11] for a discussion at zero temperature 
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and in order for the the third term in (28) to be nonsingular at the horizon, we will also 
need that 

y'iPc) = ^ (35) 

pc 

i.e. the brane should be perpendicular to the horizon. The precise of value of pc is determined 
by the boundary condition (25). 

For regularity we will also require the gauge potential vanish at the horizon^^, i.e. a{pc) = 
0. From equation (21), we can express the chemical potential in terms of density e as 

.,^A,,r,,v^, (36) 



The main technical task of the paper is to determine the behavior of fig in the limit of small e, 
which will be carried out in next section. Once the expansion of pq in terms of e is known, 
we will be able to determine the order of the phase transition and other thermodynamic 
properties of the system. This will be carried out in section IV. 



III. SMALL DENSITY EXPANSION OF THE CHEMICAL POTENTIAL 

In this section we study the behavior of (36) in the small e limit. Expanding (36) in small 
e is somewhat complicated since the solution y{p) depends nontrivially on e. The expansion 
of y{p) in terms of small e is also involved; one cannot treat the e-dependent term in (28) 
as a small perturbation of (29) uniformly for all values of p since the term becomes of 0(1) 
for sufficiently small p (when p^q^ ~ 0(e^)). 

The main result of the section is equation (60). Readers who are not interested in the 
detailed derivation can go directly to (60) and subsequent discussions. 



A. Expansion of the solution 

Since the perturbation in e is not uniform in p one needs to divide the p axis into different 
regions, treating the perturbations in each region separately, and matching them together 
at the end. For our purpose, it turns out enough to split the p-axis into two regions: 

• Inner p = e^a for ac < a < A 

Otherwise the one-form Aodt will be singular at the horizon since the norm for -Mr vanishes there. 
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• Outer Pa < p < oo 

where 



Pc = e2(Jc, pA = e2A. (37) 



The reason for the choice of scahng in the inner region can be seen by letting p = e^cr 
in (28) and one finds that for e = ^ a nontrivial scahng hmit exists, which results differential 
equation in (44) below. Also see Appendix B for a discussion of the T = solution where 
this scaling is evident. 

We will consider the limit 

e — > 0, cr = finite, A — > oo, pA — ^ . (38) 

We expand the solution in the inner and outer regions as^^ 

yx{(y) = YQ{a) + eFi(a) + e^Y2{a) + ■ • ■ (39) 

yo{p) = yoip) + mip) + ^%2ip) + ■■■ ■ (40) 

Plugging (39) and (40) respectively into (28) and expanding the equation in e, we obtain a 
series of differential equations for various functions in (39) and (40). These equations are 
rather complicated and not exactly solvable. We will work out the behavior of Yo{a), Yi{a) 
in the large a limit and yo{p)iyi{p) in the small p limit and match them in the overlapping 
region. Fortunately it turns out this is enough to find the leading expansion of p^ in the 
small e limit analytically. 

Let us first examine the outer region, yo is simply the solution to (29) which describes 
the embedding at zero density. It satisfies the boundary condition (25) at the boundary and 
near p = has the following expansion 

yoiP) = 1 + - 3(1 _ ^8)3 P + Oip') ■ (41) 

yi satisfies a homogeneous linear equation obtained by linearizing (29). For small p it has 
an expansion, 

yi = b^ip-^ + bo\ogp + h + 0{p) (42) 



While it is not a priori obvious the expansion parameter should be e and not some other power, the 
expansions below do yield consistent results. 
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where b-i and bi are integration constants and 



r^-i . (43) 



(1_ ^4)2(1 + ^4)2 

Note that the first two terms diverge in the p — hmit. The coefficient 6„i of the quadrat- 
ically divergent term will be fixed later by matching with the inner solution, hi should be 
determined by requiring yi{p oo) since yo already satisfies the required boundary 
condition (25) there. 

Now we look at the inner region. lo(c") satisfies a second order non-linear differential 
equation (primes below denotes derivative w.r.t. cr) 



^ + 3^n-^j .y„ + ^^^^pj-^ + ^^^^ = (44) 

From (34) and (35) the boundary condition at the horizon is 

ro(cTo) = r/, n'(^^o) = - (45) 

where (Tq is an integration constant, ctq is the zero-th order term in the expansion in e of 
(Tc = Pc/e^^^ (where again it is a non trivial fact that this expansion begins at 0{1) in e) 

o-c = o-Q + e 0-1 H (46) 

and can be determined by the boundary condition at infinity, 

Fo(o-) ^1 as a ^ oo . (47) 

The above condition is fixed by comparing with the leading term in (41). At large a, Yq has 
the expansion 

2(l+r/4)3/2^ +4(_l+^4)(l + ^4)4^ 

77^(27-13177^+153/7^-10577^2) . „^ „ ^ ^^^^ 

Note that the coefficient before cr~® is an integration constant which can be fixed by (45). 
We will not need its value below. 

Y'i((t) satisfies a second order linear differential equation with coefficients that depend on 
lo(c")- It is not homogeneous and rather complicated. We will not write it here. For large 
a we find the expansion 

Yi{a) = a_ia^ + aologa + ai + 0{a'^) (49) 
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where oi is an integration constant (the other integration constant presumably appears at 
higher order in the expansion) and 

a_i = r, ao = T . (50) 

Note the first two terms of (49) are divergent as a ^ oo. The constant ai can be deter- 
mined by matching with the solution of the outer region. The other boundary condition is 
determined by the regularity condition (34-35) (expanded to this order) at the horizon, this 
gives only one boundary condition on Yi and its derivative because the first order correction 
0"! in (46) is a free parameter. This parameter ai would then be determined by the resulting 
solution. 

We will now determine various integration constants above by matching the solutions in 
the inner and outer regions around pA. Note given the relation p = e^a, the e expansion 
in the outer region is reshuffled compared to that in the inner region. Comparing the first 
term in (42) to the second term of (48) we find that 

1 14r/^ 



I 2 



where we have used (43) in obtaining the second equation. From (51) and (50) we see 
that the logarithmic term in (42) precisely match with that in (49); this is a nontrivial 
self-consistency check of our expansion. It now remains to match the constant terms in (42) 
and (49). Given the difference in the argument of logarithmic term in (42) and (49), we 
thus find that 

1 7r7^ 

oi = 6i + -feologe = &i + — -yloge. (52) 

2 (1 _ ^4)2(1 ^_ ^4)2 

As remarked below (43), bi is determined by the boundary condition yi(oo) = and since 
the equation for yi is independent of e, so is fei.^'^ Thus we see from (52) that the integration 
constant ai for inner solution Yi now contains a loge piece!^^ This will be important in our 
determination of the leading order behavior for fig below. 

As a check on our above results we can make sure that the exact (numerical) solution 
at small densities matches well onto our expansions in the two different regions. This is 
demonstrated in Fig. (5). The agreement is very precise suggesting the two regions we have 



Determining bi requires solving the two-point boundary value problem for yi. 
This in turn implies that <Ji in (46) also contains log e. 
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used are sufficient. Higher order corrections in the two regions would appear as divergences 
toward the exact solution. 

y 



1.05 
/ 

y 


/ 








0.95 












0.90 












0.85 












0.80 






-10 


-5 




5 



Log(p) 



FIG. 5: The exact embeddings for e = 10^ , 10" , 10" shown in green. To clearly show the two 
regions we have plotted y{p) as a function of log(p). These are compared to the zeroth order inner 
(solid) and outer (dashed) solutions. The inner solution YQ(a) is a function of o" = pe~^^'^ so for the 
various fixed densities the curve is simply shifted along the log(/>) axis. The agreement between 
the three curves on the overlapping regions is very precise. 



B. Expansion of the chemical potential 

We now look at the small e expansion of (36). We will give the main steps and final 
results, leaving further details to Appendix B. 

The fig integral (36) can be split into integrals over the inner and outer regions 

N= (■■■)+/ {■■■) = l^T + f^o . (53) 

J Pc J PA 

fij and fia can now be expanded in terms of e using (39) and (40). The outer region 
contribution fio starts with order e, i.e. 

fio = efi^i^ + 0(e^) (54) 

and for ^u^'' only yo is needed. Using the small p expansion of i/q (41) we find that the 
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integral for /i^"* contains a quadratic and a logarithmic divergent term in the limit pa ^ 0, 

where Kq denotes the part which remains finite in the limit Pa ~* Kq, whose explicit 
expression is given in Appendix B, depends on the full solution yo{p) and can only be 
evaluated numerically. 

The inner region contribution also contains an 0(1) piece 

/ix = /if + e4'^ + ■ ■ • (56) 

with (recall that a = e^^p) 

(■■■)-/(...) (57) 



where in the second equality we have separated the expression into two pieces by splitting the 
integral. Now changing the integration variable of the first integral into Yq and comparing 
it with (26), we find that it is precisely m''^\ The second term can be evaluated in power 
series of 1 /A by using the expression (48) for Yq at large a, leading to 



,(0) Lq (1 - 

To evaluate jjSp one also needs Yi. Using the expansion (49) of Yq and Yi for large a, we 
find that the integral for contains a logarithmic divergence (in the A — oo limit) from 
the upper end of the integral and can be written as 

(1) Lo /2r/4(3-r/^ + 3r/8) , ,., Tr/^loge , , , , .\ ,,,, 



27ra' V(l + (1 -77^) (1 + 7^4)4 ^ 

Note that in the finite part we have isolated a piece proportional to log e which comes from 
the integration constant ai of Yi through equation (52) (see Appendix B for details). Kj is 
finite and independent of e. Its explicit expression is given in B and can only be evaluated 
numerically. 

Now adding (55) and (58), (59) together we find that the second term in (58) precisely 
cancels the quadratic divergence in (55) and the logarithmic divergence in (59) precisely 
cancels with that in (55), leaving a finite piece proportional to eloge. We thus find that 

^ = m^p - B{T)e log e + A{T)t + 0{e^) (60) 
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where 

B{T) = Bi{T) + B2{T) 
and 



Lo r/^ (3 - 8r/^ + 3r/S) 



2Txa' (1 -ry4)(l + r7Y ' 



B2{T) 



A{T) 



Ln 



2na 



- {Ko + Kj) (61) 



7rf 



(62) 



27ra' (1 -774)(1 +r/4)4' ^ 2™' (1 - r/4)(i + ^4)4 

Note that in writing (61) and (62) we have emphasized that the coefficient of eloge receives 
contributions from two different sources; Bi from cancellation of the logarithmic divergence 
and B2 from the finite piece in (59). The expressions for Kq and Kj are given in Appendix B, 
which can only be evaluated numerically. We also note that in the e ^ limit it is the e log e 
which dominates and whose coefficient we have determined exactly. Note that Lq, which 
gives of the location of the tip of the D7-branes in the absence of a charged density (13), is 
also temperature-dependent. 

x(e) 



0.98 - 




log(e) 



FIG. 6: Plots of the exponent x(e) as defined in (63), for various values of the temperature. At 
zero temperature the exponent is simply 1 and all these curves should eventually limit to 1 at small 
enough densities. The slow running of this exponent is a consequence of the loge behavior of the 
chemical potential. The dashed line is a (1 parameter) fit to the numerical curve using (64). 

Again we can use exact (numerical) solutions for small densities to check our analysis. 
The loge dependence is not easy to see, so we must push our analysis to very small densities 
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with a wide range of densities. Then the best way to extract the behavior is to plot the 
function x(e) defined by, 

xie) = ^M^.-^^^) (63) 
ct log e 

^ A(T)-i^(T)(loge + l) 

~ AiT)-BiT)loge ^'^^ 

where for the last equality we have used our derived small e expression (60). We can fit this 
later form to the numerical result and find the ratio A{T)/ B{T), after which we can use the 
overall normalization of fiq — m^^^ to fix A{T) and B{T) separately. This agreement in the 
form of the solution is shown in Fig. 6. 

C. Behavior of B{T) and A{T) 




FIG. 7: The quantities in (60) are plotted as a function of T/T^. The dots represent the values 
of A{T) and B{T) obtained through fitting the numerical results to the form of the small density 
expansion. The actual curves come from the results of section III (which also require numerics to 
calculate A[T).) The consistency is gratifying. 

Let us now examine the behavior of B{T) and A{T) as we vary T. In the T ^ (i.e. 
?7 — i> 0) limit, from (61) 

m,(T) = m„ A(T)=m,y, 5(T) = (65) 
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which reduces to the correct zero temperature behavior (A4) in Appendix A. In the above 
equation k is a number defined by the integral k = dx/{x^ + 1). For general temperature 
T, the behavior of A(T), B{T) and m^"^^ are plotted in Fig. 7. Note that in obtaining Fig. 7 
we have converted the dependence on rj into T and also the temperature dependence of Lq 
using the relations demonstrated in Fig. 3. 

Note that B{T) is positive at small temperature and becomes negative for 

r/>r/„ r/, = (^(4- v^))^ i.e. T > T„ T, ^ .982Trf, fi(T,) = (66) 

Also note that A{T) becomes negative at some temperature higher than T^. Furthermore it 
can be shown that A{T) diverges at a temperature (TmLOOST^), beyond which Minkowski 
embedding no longer exists; clearly our analysis does not apply since there is no zero density 
solution about which to perturb in e. 



IV. THERMODYNAMICS 



From our small e analysis we can now extract some important aspects of the thermody- 
namics of this gauge theory. We will mostly work in the grand canonical ensemble, where 
we fix T, /i^ and use the pressure -P(/ig, T) as the appropriate thermodynamic potential. 



A. A third order phase transition for < T < Tc 

As noted above B{T) > for < T < Tc. In this case (60) may be inverted to find 

eiyiiq) which is then a single- valued function of Hq for Hq > In other words there is 

only a single black hole embedding for a given fixed fiq > Thus for B{T) > we 

(T) (T) 

conclude there is a continuous phase transition at /Xc = rriq ■ For /ig < rriq we have a 
single Minkowski embedding with e = 0. 

Examining derivatives of the pressure on either side of rriq one finds for /i^ < ruq , 

P(T, /i) = Po(T), nq = -^ = 0, ^ = 0, for all n > 1 (67) 



and for fiq > m}[\ we find 



d^q d^f' 



g^^o^fg^y, (68) 



(9/i2 \de J A-B-B\oge 
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and 

a/i3 de^ \dej e{A-B-B\ogey 
In the limit e ^ 0, — > and ^ oo. This imphes that there is a third order 
phase transition since the third derivative of the pressure is discontinuous across the phase 
boundary. Also since > 0, the system is thermodynamically stable. 

Note that exactly for B(T) = the transition becomes second order, since (68) is then 
discontinuous across the phase boundary. This happens at two places, the first is the zero 
temperature critical point which was studied in [11]. The second is a new critical point at 
(T, /i) = (Tc, /ic) where Tc is given in (66) and fic = ~ ~ .0418mg. We will discuss the 
phase structure around this point in the next subsection. 

In Fig. (1) we display the phase diagram in the fiq — T plane, where this 3rd order phase 
transition is the predominant feature. 



B. First order phase transition for T > Tc 

Beyond the critical point T > Tc we find that B{T) becomes negative. From (68) this 
implies that becomes negative in the limit of small e indicating a thermodynamic insta- 
bility. 

To see what happens beyond Tc, let us plot fig at fixed T as a function of e. As indicated 

(T) 

in the left plot of Fig. 8 the curve drops below jjq = rriq . Since we expect the curve will 
go up again for sufficiently large density, there must be at least one minimum somewhere, 

(T) 

which we call will /Xmin- This minimum satisfies /imin < ■ That is there will be multiple 
black hole embeddings for a fixed fig > /imin and which appear as a function of fig before one 

(T) 

reaches jj,g = rrig . This allows for the possibility of a first order phase transition from the 
Minkowski embedding with e = to a black hole embedding with e = €d{T) ^ for some 
chemical potential fig = fid(T) between /imin < f^d < f^g ■ In this way one circumvents the 
thermodynamically unstable solution. 

Although this discussion was quite general for T > Tc, we can make it more explicit 
by considering temperatures |T — Td <C Tc, then the minimum and the other black hole 
embedding still lie at small enough values of e that (60) applies with only small corrections. 
The left plot of Fig. 8 gives the behavior of the chemical potential as a function of quark 
density near and the right plot of Fig. (8) demonstrates the swallow tail behavior of the 
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0.0004 r 




FIG. 8: (left) The quark chemical potential as a function of quark density, for three temperatures 
T/Td = .980, .995, .996. The first (long dashes) is below Tc and has no minimum. For the other 
two we have indicated the chemical potential at which the transition occurs, (right) The grand 
potential density (related to the pressure hy P = —g) as a function of chemical potential for the 
same temperatures as the left figure. 

pressure for this first order transition. Using (60), the first order transition point fid and the 
discontinuity in the density can be expressed in terms of the functions A{T) and B(T) as 

q(T) ^ exp /i.(T) - mf ) ^ ii?e.(T) (70) 

and very close to Tc we can expand A{T) and B(T) 

A = Ac + 0{T -Tc), B = -B'^{T -Tc) + 0{T -Tc)'^ (71) 

to find 

6,(T) ^ exp = exp(-.12T,/(T - TJ) (72) 

In Fig. (9) we plot and fid as functions of T near Tc using (70). 

To conclude this section we note that our analysis near Tc relies on the small density 
behavior (60). An implicit assumption in our discussion is that various curves in the left 
plot of Fig. 8 does not turn back down at larger densities. While our small e analysis itself 
does not rule out this possibility, we believe it does not happen for temperature sufficiently 
close to Tc from our own numerical study of embedding solutions and from the finite density 
analysis of [16].^^ Our analysis of the first order phase transition cannot be trusted for 

At some higher temperature, this could happen. 
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FIG. 9: {left) The first order phase transition line (solid) close to the critical point, in the {^q 



m. 



(T) 



) — T plane. Also shown (dashed) is the region where multiple embeddings are available at 



fixed fiq. (right) The behavior of €d{T) above the critical temperature The shaded region represents 
the onset of a thermodynamic instability. In the canonical ensemble this region is circumvent as 
usual by the Maxwell construction. 

temperatures too high above - since then the density at the first order phase transition 
will lie out of the range of applicability of (60). Note that (60) is likely to break down much 
earlier than simply when e ~ 1. The reason for this is that at some temperature multiple 
embeddings for a fixed density should appear, something which our analysis has not allowed 
for. This is clearly the case (as discussed earlier) for zero density, where there exists a 
temperature T^h < above which multiple embeddings are allowed (the original Minkowski 
embedding and two new black hole embeddings). One expects this to be the case for 
small densities. Indeed these new embeddings will be continuous deformations of black hole 
embeddings away from zero density, so our starting point (Minkowski embeddings) is not 
good for seeing these multiple embeddings. As we will discuss later the reason for this break 
down is closely related to how this line of first order phase transition connects to the zero 
density dissociation transition. 



Note that T^h ~ 0.9975rd is far enough away from Tc = .982Td that one can trust that multiple embeddings 
do not appear then. 
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V. CONCLUSIONS AND DISCUSSIONS 



In this paper we showed that in the plane of temperature v.s. baryon chemical potential 
there is a critical line of third order phase transition which ends at a tricritical point after 
which the transition becomes first order. The critical behavior at the critical line is given 
by (60) which contains an intriguing logarithmic behavior. It would be interesting to have 
a microscopic understanding of this behavior, and more generally the structure of the whole 
phase diagram. In particular, it would be interesting to see whether the logarithmic behavior 
is related to the spiral behavior observed in [7, 20] at zero chemical potential. Below we 
discuss in more detail two open issues of our investigation. 

A. Connection to the dissociation transition 

To complete the phase diagram, we need to address whether the first order phase tran- 
sition slightly above discussed in sec. IV B connects to the first order phase transition 
along the vertical axis at T = and fig = 0, and if yes, how. Note that such a conclusion 
is suggested by numerical work in [10] and is consistent with the numerical work in [12]. 

For convenience of discussion let us first briefly review some important aspects of the 
transition at fig = 0. At a low temperature there is only Minkowski embedding. As one 
increases T to a temperature Tbh, two new black hole embeddings appear and when one 
further increases the temperature to > T^h = 0.997Trf, a first order phase transition 
occurs, above which one of the black hole embedding thermodynamically dominates over 
the Minkowski embedding. For all the embedding solutions, the baryon density is zero. 
Thus the discontinuity ed(T) at the phase transition is trivially zero. 

We now return to the question at hand, which unfortunately our analysis for small density 
cannot directly answer. To see this if we were to extrapolate our analysis of the 1st order 
phase transition line (the solid line in the left figure of Fig. 9), we find that the line fid{T) 
crosses the fig = axis at some temperature. This temperature is slightly higher than 
Td and hence would seem to give an estimate of itself, subject to e corrections. This 
small difference from might not seem like a problem, and one might guess we have a 
nice description of the physics of the zero density dissociation transition. However it is 
qualitatively wrong, since we also find €d{T) grows monotonically along this line. This is 
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qualitatively wrong because on the = 0, the phase transition should occur at zero 

density. This demonstrates that (60) must have broken down quite drastically somewhere 
before HdiT) = is reached. 





FIG. 10: The qualitative behavior of //g(e) at different temperatures in order for the 1st order 
phase transition near discussed in last section to connect to the transition at /ig = 0, e = 0. Left 
plot: for a temperature T slightly above T^,. Middle plot: for a temperature T above Tc and below 
Tbh. Right plot: for T >Tbh. 

We will now consider the simplest possibility, i.e. assuming that 1st order transition 
line from Tc does smoothly connect with the transition at zero fiq. For this to happen the 
transition density erf(T) should first increase as we increase T from (consistent with the 
behavior derived from (60)) and then decrease to zero as T^ is approached. In turn this 
implies that the minimum of the curve fiq{e) should approach the e = axis as the tem- 
perature is increased to T = Th^. In Fig. 10 we plot qualitatively, at different temperatures, 
what is required of the curve /ig(e) for this to happen. In particular, the minimum of the 
curve will hit the origin e = /ig = 0, exactly at the temperature T = Tbh, where two black 
hole embeddings appear at zero density. For temperatures above T > Ti,h moving towards 
T = Tfi the curve /ig(e) is similar to that at T = Tb/,. However now the point at which the 
1st order phase transition occurs (which we have called {fid,^d)) should move towards the 
origin. Note that plots of form of Fig. 10 were found numerically in [9]. 

An important feature of the last two plots Fig. 10 is that when temperature is sufficiently 
high, /i can be multi-valued for a fixed e. For T > T^h., this is consistent with our expec- 
tation that there should be three black hole embeddings at small e, two from small density 
perturbations of the black hole embeddings at zero density one from the small perturbation 
of the Minkowski embedding which we discussed earlier in this paper. We have done some 
preliminary study of the behavior of the fi — e plot near the origin for T ~ T^h from per- 
turbing two black hole embeddings at zero density and confirmed the qualitative behavior 
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presented in the last two plots. But the extrapolation to the branch we studied in this paper 
requires studying finite density solution and that part of the curve is pure speculation at 
the moment. We leave the full exploration to future work. 

Since, for high enough temperatures, there are multiple embeddings with different chem- 
ical potentials and the same density, if we examine the thermodynamics in the canonical 
ensemble, there will be an interesting set of phase transitions unrelated to the one identified 
above that can occur between embeddings with different chemical potentials, but the same 
density. These phase transitions are quite strange, and as discussed in [16] occur between 
solutions which are potential thermodynamically unstable. These transitions should involve 
embeddings with chemical potentials higher than that at which the 1st order phase transi- 
tion in the grand canonical ensemble occurs. Hence such transitions will be hidden from the 
point of view of the grand canonical ensemble. Since thermodynamics should be consistent 
in the two different ensembles one would expect that a proper consideration of the canoni- 
cal ensemble using for example, the Maxwell construction of inhomogeneous phases, should 
remove these differences. 



B. Transition at finite A 
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FIG. 11: Possible phase diagram at finite A. The lower region consists of a Boltzmann gas of 
quarks and anti-quarks. The third order phase transition of Fig. 1 is potentially smoothed to a 
cross over. 



Another important question is what could happen to the phase diagram at finite 't Hooft 
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coupling A. At finite A, as pointed out in [21] and commented upon earlier in [12], for any 
/ig 7^ at any finite temperature, the baryon charge density is in fact nonzero, given by a 
dilute Boltzmann gas of quarks. Thus in the phase diagram Fig. 1, the region below the 
critical line also has a finite density. This is of course expected for a deconfined plasma. In 
particular the value of quark density at the transition line should be nonzero and continuous 
at finite A. This immediately raises the possibility whether the transition is smoothed into 
a crossover at any finite A. To have a definite answer to this question requires summing 
over the instanton contributions in (32) for all n and then taking the limit /i^ — m"q^ in the 
resulting expression. This is certainly beyond the scope of the current paper. In Fig 11 we 
plot qualitatively the structure of the phase diagram if the transition is smoothed out. 

Here we mention another indirect indication that the transition might be smoothed out. 
In [21] we have argued that at any finite A, in the region below the transition line in the 
phase diagram, the mesons have a width proportional to the sum of quark and anti-quark 
densities, 

r = ^f^l^(^ = 0)rK + n_). (73) 

L q 

where n± oc exp(— (mg"^-* ± jj,q)/T) and are exponentially small in \/X as argued in sec. II C. 
iIj{6 = 0) is the meson wave function at the tip of the brane. We have also performed [22] a 
calculation of the meson width on the critical line from above and found at small densities 
exactly the same answer as (73).^^ Note that one of the key signatures of the phase transition 
at zero density (in fact potentially the defining signature) is the spectrum of mesons (quark 
and anti-quark bound states) [19, 20]. Here we find the meson widths and thus the spectrum 
are continuous across the critical line indicating that the transition might be smoothed out. 
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APPENDIX A: THE ZERO TEMPERATURE SOLUTION AT SMALL DENSI- 
TIES 



One way to motivate the scahng of the inner region p = cre^/^ is to look at the zero 
temperature solution in detail, where the exact solution is known. 
This solution is: 

y{p) = \ , a'ip) = -y'ip) (Al) 

where c is an integration constant related to the expectation of the mass operator dual to 
the field y on the gravity side. The boundary conditions are y{0) = a(0) = 0. 

We can fix c in terms of the quark mass, or equivalently in terms of L. After rescaling 
as (23) hj Lq = L the condition y{oo) = 1 fixes the quark mass. Integrating (Al) one finds 
for c, 

«:3c3 = (^2 _ c2) 

where k is a number defined by the integral, k = dx/{x^ + 1). For small densities the 
condensate c approaches e as, 

c^e- {n/Lfe^ , - ^ e^{K/Lf (A2) 



Note that by small densities we really mean e ^ or equivalently nq/(NcVX) <^ m^. 

From (Al) it is clear the cross over scale for y{p) between a flat Minkowski embedding 
and the curved solution going to the (AdS) horizon is ~ — or p ~ e^^^ at small 
densities. This is the scaling that we set out to demonstrate. 

Rescaling a = p/e^^^ we find the zero temperature inner solution becomes. 

In this limit the chemical potential to order e receives contributions only from inner region 
with the result 

r. / 1 _ _ \ 

(A4) 
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APPENDIX B: DETAILED ANALYSIS OF THE CHEMICAL POTENTIAL 



In this appendix we give some details in the small e-expansion of fig, equation (36), which 
we copy here for convenience 

Various functions in (Bl) were introduced in (24). We split the above integral into those 
over the inner and outer regions 

N= / (■■■)+ / (■■■)=/ix + /io (B2) 

J pc J PA 

which can then be expanded in terms of e using the expansions (39), (40) of y{p) in inner 
and outer regions respectively. 

Let us first look at fio- We find that it starts with order e 

/io = e/ig) + 0(e2) (B3) 

with 



and = + Voip)- Using the expansion (41) of yo at small p, we find that the integral 

contains a quadratic and a logarithmic divergent term in the limit p\ ^ 0, given by 

2ira' ,1, 1 - rf 2r,*{3 - rf + ?.r,») , 

-L^"" = 2(l+.,ypi - (1 + ,Y + *° 

where i^o denotes the finite part of the integral. It can be defined more precisely by 



Ko = lim 

PA-*0 



(1 - r^Vu') y/IT^ (1 - v') 2rj\3 - + 3^') 

PA ^ il + vVu^fp' ^a + V^Pl (-l + ^^)(l + r/T 



(B6) 

Once yo{p) is known numerically, Ko can be calculated numerically. 
We now look at px which can be expanded as 

Pi = pf + e/i?^ + ■ ■ • (B7) 

The first term in the above has been discussed explicitly in the main text. Here we give 
some details for computing /ij"*, which can be written as (derivatives are w.r.t. a) 



2W m f'' (Yft-V^) 

-4\l/2 
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where 

The total derivative term in (B8) is nonzero only at the upper end and gives 
total = , „ A- + (ifjfe + '>"(3 + .>-)-lVlog(A) 

where we have used the expansions (48) and (49). The integration over M in (B8) is also 
divergent in the limit A ^ oo. The divergent terms can be extracted using the expansion (48) 
and one finds that 

^ rf^ M = -—Jf-—A^ - 2r/^(3- 8r^^ + 3r/2 

(1+7/4)3/2 (_1+ ^4)^1 + ^4)4 SI i K J 

where Mf denotes the finite part of the integral and can be defined more precisely as 

M, = hm ([\aM+ / + f^^' ^^ll^l ^^^(A)) (B12) 

Adding (BIO) and (Bll) together we find that the quadratic divergences in A cancel and 
27ia' (I) 

V(3-r/^ + 3r/S) Ir^Hoge 1 - r/^ 7/^^(3 + 7/^) 

(l_^4)(l+^4)4 ^+(l_^4)(l+^4)4+^/^^ 1 (l_^4)(l + ^4)4+ ^ 

(B13) 

where we have used (52) and hi was introduced in (42). Equation (B13) gives (55) with the 
definition 

(B14) 

Kj depends on the full solution and can be computed numerically. 
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